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(57) An adaptive image tone mapping curve (170) 
based on perceptuai preference guidelines is generated 
as a sigmoidai function, in which the sigmoidai function 
parameters (150) (slope and shift) are deternnined by 
original image statistics (149). Tone curves (170) gen- 
erated for different Images each have a smooth sigmoi- 
dai shape, so that the tone mapping process does not 
change the Image histogram shape drastically. The sig- 
moidai function has the form: 



t{x) = 



100 



where a is the slope parameter (1 50a) and |3 is the shift 
parameter (150b). The input value x in the sigmoidai 
function varies In the range [0, 100], because the tone 
curve (170) is generated on an L* scale, which has val- 
ues from 0 to 1 00. The sigmoidai tone curve (1 70) cal- 
culation can be implemented efficiently using simple 
arithmetic operations by pre-calculating and storing var- 
ious factors used In the calculation of a and ^ and by 
pre-generating a pair of fixed tone curves (158) with two 
extreme slopes and Interpolating between the curves 
(158). 



1 + exp(-a(x/1 00 - P)) ' 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001] The present invention relates generally to digital innage tone mapping algorithms, and specifically to digital 
image tone mapping algorithms based on perceptual preference guidelines. 

^0 Description off Related Art 

[0002] Digital image tone mapping referstotheprocessof mapping captured scene luminance levels to the luminance 
or densrty levels of an output device (display or printer). Digital image tone mapping is necessary due to the fact that 
scene luminance ranges very seldom match the luminance range of the output device. Digital image tone mapping 
t5 algorithms can be implemented within any digital device to produce visually pleasing output images. For example, such 
devices car include digital cameras, scanners, digital camcorders, printers capable of printing digital images and digital 
televisions. In addition, digital tone mapping algorithms are frequently used in image enhancement software applica- 
tions. 

[0003] Since the human eye adapts to different light levels in different ways, perceptual factors must be taken into 
20 account When attempting to generate a rendered image that iooks "right" with respect to the original scene. Depending 
on the application, image tone mapping typically has one of three goals; (1) appearance match; (2) subject ve prefer- 
ence; or (3) Information preservation. The goal of appearance matching strives to make the rendered image as per- 
ceptually similar as possible to the original scene. This is usually an implicit goal In consumer imaging and image 
synthesis applications. Taking into account subjective preferences allows the image to look as pleasing as possible to 
25 the viewer. This is usually desirable in consumer imaging and commercial photography. If the goal is information pres- 
ervation, the algorithm seeks to preserve or enhance the details at all regions and all luminance levels of an image. 
This is most often requested in medical imaging, satellite imaging, and archiving. 

[0004] Many existing tone mapping algorithms focus on achieving an appearance match between the original image 
and the rendered Image. There are at least two perceptual factors typically considered In such algorithms: (1) global 
^0 luminance adaptation; and (2) local luminance adaptation. The overall global luminance level of the scene Influences 
the adaptation state of the eye. Two aspects of such global luminance adaptation have significant impact on tone 
mapping, brightness and spatial contrast. 

[0005] First, the brightness function Is different at different scene luminance levels. Perceived brightness, which 
corresponds to the viewer's perceived strength of light, is roughly a power function of physical luminance (Steven's 

35 law). The exponent for such a power function is larger when the overall luminance level is higher. When rendering an 
image on a lower luminance device, the exponent must be adjusted to accommodate such differences. In addition to 
the brightness function change, the spatial contrast sensitivity of the eye also changes as it adapts to different ambient 
luminance levels. When the ambience is bright, the eye perceives the high spatial frequency components (details) of 
an image better than when the ambience is dark, i.e., the visual acuity of the eye improves with better ambient lighting. 

40 Also, the contrast threshold. I.e., the minimum contrast needed to detect components of the image, decreases with 
Increased luminance level. To render a bright Image onto a lower luminance device, the luminance contrast of the 
details In the Image can be enhanced to account for these effects. 

[0006] Two different tone mapping algorithms developed by Jack Holm and Tumblin & Rushmeier, respectively, ac- 
count for the brightness function change by adjusting the curvature of the tone curves based on the absolute luminance 
level of the scene, Each of these tone mapping methods Is described separately In Holm, J., "Photographic Tone and 
Colour Reproduction Goals," CIE Expert Symposium on Colour Standards for Image Technology, pp. 51-56 (1996); 
and Tumblin, J. and Rushmeier, H., "Tone Reproduction for Realistic Images," IEEE Computer Graphics and Applica- 
tions, 13(6):42-43 (1993), both of which are hereby incorporated by reference. These algorithms have the benefit of 
creating the proper overall sensation of brightness or darkness corresponding to the original image, which is desirable 

50 in high-end digital imaging. However, both require accurate information about the absolute luminance level of the 
original image. In a digital camera, it is possible to estimate absolute luminance levels of image pixels from the raw 
pixel values and the camera's capture settings such as aperture, exposure time, lens properties, etc. However, In low 
cost cameras, such calculations are often not available due to the added complexity and cost. Therefore, a global 
luminance adaptation solution that focuses on preference is more practical. 

55 [0007] In tone mapping algorithms based on preference, the goal is to achieve a certain set of image properties liked 
by viewers. The widely used histogram equalization method can be categorized as such an algorithm. The histogram 
method is based on the observation that most "good" images have a luminance histogram that fully occupies the output 
dynamic range. The algorithm adjusts image gray levels to move the histogram shape toward a flat, Gaussian, or some 
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other predetermined shape. Of course, how well such a method works depends on whether the assumption Is true that 
every "good" Image has the same histogram. The method does well on Images that have a symmetric and well-dls- 
tributed histogram, but makes Images look unnatural when there are large areas of dark or light background In the 
image, which bias the histogram toward one side. 

5 [0008] A modified histogram equalization method developed by Larson, etal. is more robustthan traditional histogram 
equalization methods. Larson's method limits the amount of gray level adjustments allowed In the tone mapping. The 
amount of gray level adjustments are limited based on luminance contrast sensitivity measurements. In addition, one 
variation of this method also accounts for the change in visual acuity under different Illumination levels. Reference is 
made to Larson, G. W., Rushmeier, H. and Piatko, C, "A Visibility IVIatctiing Tone Reproduction Operator for High 

^0 Dynamic Range Scenes," IEEE Transactions on Visualization and Computer Graphics, 3(4):291-306 (1997), which is 
incorporated by reference. However, the computation Is Iterative, and thus the implementation is costly and slow. In 
addition, the modified histogram equalization method also requires accurate absolute luminance level information. 
Thus, although this modified histogram equalization method creates a more accurate appearance match, it does so 
at the cost of higher computational complexity. 

15 [0009] The local luminance adaptation perceptual factor considers the fact that the eye looks at an image by scanning 
around. The eye can rapidly adapt to the luminance level of small regions in the orig nal scene to enable regions in 
the shadows and in the highlights to be clearly visible to the eye. In the rendered image, both the dynamic range and 
the adaptation environment are different. Therefore, to fully Imitate the eye's adaptation process, the luminance levels 
of an image are adjusted according to its local luminance levels. 

20 [0010] Various local tone mapping algorithms, such as Tumblin's detail-preserving contrast reduction method (Tum- 
blin, J. and Turk, G., "LCIS: A Boundary Hierarchy for Detail -Presen/ing Contrast Reduction," Computer Graphics 
Proceedings, SIGGRAPH 99, pp. 83-90, Los Angeles, CA, USA (1 999), which Is incorporated by reference), and various 
algorithms based on the retlnex theory have attempted to Imitate the local luminance adaptation process. Reference 
Is madeto Jobson, Rahman, Z, and Woodell, G., "A Multlscale Retlnex for Bridging the Gap Between Color Images 

25 and Human Observation of Scenes," IEEE Transactions on Image Processing, 6(7):965-976 (1997); and Rahman, Z., 
Jobson, D. and Woodell, G., "Multi -Scale Retlnex for Color Image Enhancement, Proceedings," International Confer- 
ence on Image Processing, volume 3, pp, 1003-1006, Lausanne, Switzerland (1996) for a discussion of the retlnex 
theory, both of which are incorporated by reference. Although these algorithms do preserve the local contrast of images, 
they are Iterative methods that also Involve the decomposition of different spatial resolution components of an Image, 

so which Is computationally costly. 

SUMMARY OF THE INVENTION 

[0011] A system and method for Implementing an adaptive digital Image tone mapping algorithm that Is based on 
35 perceptual preference guidelines is provided. The tone curve is generated as a sigmoldal function having sigmoldal 
function parameters (slope and shift) that are determined by original image statistics. Since the parameters are based 
on original image statistics, the algorithm is adaptive. Tone curves generated for different images each have a smooth 
sigmoldal shape, so that the tone mapping process does not change the image histogram shape drastically. Therefore, 
the algorithm is both robust and conservative (e.g., while Improving the appearance for most Images, the algorithm 
40 does not make any Images look worse). The sigmoldal function has the form: 

1 +exp(-a(x/100-p))' 

where a Is the slope parameter and |3 Is the shift parameter. The Input value x In the sigmoldal function varies In the 
range [0, 1 00], because the tone curve Is generated on an /. * scale, which has values from 0 to 100. The sigmoidal 
tone curve calculation can be implemented efficiently using simple arithmetic operations by pre-calculating and storing 
various factors used in the calculation of a and p and by pre-generating a pair of fixed tone curves with two extreme 
slopes and interpolating between the curves. 

50 

BRIEF DESCRIPTION OF THE DRAWINGS 

[001 2] The disclosed invention will be described with reference to the accompanying drawings, which show important 
sample embodiments of the invention and which are incorporated in the specification hereof by reference, wherein: 
55 [0013] FIGURE 1 is a block diagram Illustrating a system utilizing a preferred embodiment of the tone mapping 
algorithm of the present invention; 

[0014] FIGURE 2 Is a flow chart Illustrating the steps Involved In the tone mapping algorithm shown In FIGURE 1 of 
the drawings; 
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[0015] FIGURES 3A-3C are blockdiagrams Illustrating alternative embodiments for implementing t^^ 
processor shown In FIGURE 1 of the drawings; 

[0016] FIGU RES 4A-4C are flow charts Illustrating the steps involved in eacli of the embodiments shown in FIGURES 

3A-3C, respectively: 

5 [0017] FIGURE 5 is af low chart illustrating thesteps for calculating the slope parameter of the tone mapping algorithm 
of FIGURE 1; 

[0018] FIGURES 6A and 6B are flow charts illustrating alternative embodiments for calculating the shift parameter 
of the tone mapping algorithm of FIGURE 1 ; 

[0019] FIGURES 7A and 7B are biocl< diagrams Illustrating alternative embodiments for generating the tone curve 
^0 of the tone mapping algorithm of the present invention; 

[0020] FIGURES 8A and 8B are flow charts Illustrating the steps Involved In the alternative embodiments shown In 
FIGURES 7A and 7B, respectively; 

[0021] FIGURE 9 is a graphical representation of a sample tone curve of the type generated by the tone mapping 

aigorithm of the present invention; 
15 [0022] FIGURE 10 is a graphical representation of two pre-generated tone curves plotted on L*axes, in accordance 
with the embodiment shown in FiGURE 7B and 8B; 

[0023] FiGURE 11 A is a graphical representation of output L'S' levels plotted as a function of the slope parameter, in 
accordance with the embodiment shown In FIGURE 7B and 88; and 

[0024] FIGURE 1 1 B is a graphical representation of gamma corrected output levels plotted as a function of the slope 

20 parameter, In accordance with the embodiment shown In FiGURE 7B and 8B. 

DETAILED DESCRIPTION OF THE PRESENTLY PREFERRED EXEMPLARY EMBODIMENTS OF THE PRESENT 
INVENTION 

25 [0025] The numerous innovative teachings of the present application will be described with particular reference to 
the presently preferred exemplary embodiments. Howeven it should be understood that these embodiments provide 
only a few examples of the many advantageous uses of the innovative teachings herein, in general, statements made 
in the specification of the present application do not necessarily delimit any of the various claimed inventions. Moreover, 
some statements may apply to some Inventive features, but not to others. 

so [0026] With reference now to FIGURE 1 , a digital Image tone mapping aigorithm 1 00 in accordance with preferred 
embodiments of the present invention is shown implemented within a digital image system 1 0. The digital image system 
1 0 can be any digital device, Including, but not ilmlted to, a digital camera, scanner, digital camcorder, printer capable 
of printing digital images or digital television. The aigorithm 100 can be included within an image processing chip of 
the digital device, such as a low-end CMOS sensor chip, a CCD sensor chip or a separate image processing chip 

35 packaged together with the sensor chip, or within a software package (driver) associated with the digital device. Alter- 
natively, the digital image system 10 can be a computer system having an Image enhancement software package 
installed thereon in order to enhance tonal characteristics of digital images. 

[0027] The tone mapping algorithm 1 00 is simple and fast in terms of memory and processing time, and uses a global 
tone mapping method, with one tone curve for ail color channels. It should be noted that the tone mapping aigorithm 

40 1 00 can be applied to both color Images and monochrome Images. The tone mapping algorithm 1 00 produces a smooth 
and consistent tonecurve^ with gentle curvature, which avoids the extreme contrast changes evident in some histogram 
equalization methods, i n addition , the tone mapping aigorithm 1 00 is ideal for video applications since the tone curve 
remains relatively stable from frame to frame by using an adaptive tone which is always sigmoidai in shape. Further- 
more, the tone mapping aigorithm 1 00 generates images that look pleasing to the user (in terms of brightness, contrast 
and color constancy) without requiring an estimation of the absolute luminance levels of the original Image. 
[0028] The tone mapping algorithm 100 is based on perceptual preference guidelines, such as the preservation of 
general histogram shape at the middle range by using a centering function on L* nth and (1 00-/i)th percentile values, 
an standard deviation of around 20 and a mean L* vaiue of around 50. it should be understood that the scale is 
simply a transformation from a linear scale based on received linear image data, It should be noted that other perceptual 

50 preference guidelines can be used based on the requirements of the user or the specific digital image system 1 0 in 
which the tone mapping aigorithm is implemented, in order to produce atone mapping algorithm 100 that takes into 
consideration perceptual preferences and that compiles with computation complexity constraints, a sigmoldal function 
tone mapping method is used, which will be described in more detail hereinbelow with respect to FIGURES 2-8. 
[0029] With reference now to FIGURE 2, which will be described in connection with FIGURE 1 , there Is illustrated 

55 the steps for mapping received digital image data representing the luminance (pixel) values of a digital image to the 
luminance range of an output device 50, using the tone mapping algorithm 1 00 of the present invention. A digital image 
source 20, such as a digital image capture device or a digital image file, provides digital Image data 25 associated with 
a digital image to a linear transformation device 30 (step 200), which transforms the digital Image data 25 Into linear 
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Image data 130 on a linear scale (step 21 0). It should be noted that the original digital Image data 25 may be on any 
scale. Alternatively, the digital Image data 25 may already be on a linear scale, which obviates the need for the linear 
transformation device 30. This linear image data 130 Is later used by a tone mapping operator 40 to map tlie linear 
image data 130 onto the luminance range of the output device 50 (step 250). 

5 [0030] Once transfonned, til e linear image data 130 is passed to an image statistics processor 120, which determines 
one or more Image statistics 149 based on the original linear Image data 130 (step 220). These Image statistics 149 
are used by tone curve parameter calculation logic 140, along with one or more predefined perceptual preferences in 
the calculation of two tone curve parameters 150, namely, the slope parameter and the shift parameter (step 230). The 
pre-defined perceptual preferences may be defined by an operator (e.g., a user or manufacturer of the digital Image 

10 system 10). The shift parameter shifts the tone curve so that under or over-exposed images are brought closer to the 
rigtit range. The slope parameter stretches or compresses the tone curve to achieve a desired histogram spread. These 
two tone curve parameters 150 are used by a tone curve generator 160 in generating a tone curve 1 70 for the image 
(step 240). The tone curve generator 1 60 produces a sigmoidal tone cun/e 1 70 using the following sigmoidal function: 

100 

K^) = 1 ; — 7—FT7^ — (Equation 1) 

^ ' 1 +e){p(-ot(x/100 - p)) ^ ^ ' 

where a is the slope parameter and |3 is the shift parameter. The input value x in the sigmoidal function varies in the 
range [0, 1 00], because the tone curve 170 is generated on an scale, which is a measure of the brightness level of 
20 the image pixels and has values from 0 to 100. Based on the CIE standard formula, the /_* scale is approximately a 
cube root transformation of the linear image data: 



25 



L* = 116(^)^ - 1 6, (Equation 2) 



where Y\s linear luminance, and is the luminance level of the whitepoint. Due to thefactthat lvalues are normalized 
by Vp values for the calculation of the absolute scale of Vand Y„ is not important, as long as both are proportional 
to the absolute luminance. In implementation, the raw or color-corrected linear RGB values may take the place of the 
50 V values. 

[0031] The L * scaie is the preferred scale in the generation of the tone cun/e 170 for several reasons. First, the L 
* scale Is perceptually more uniform than a linear scale. Second, most natural Images tend to have more symmetric 
histograms on an L * scale than on a linear scale, which makes the adjustment of the histogram spread easier (and 
perceptually more robust) on the L * scale. However, it should be noted that a linear scale may be used by the tone 
35 curve generator 160 instead of the L * scale described herein. Also, it should be understood that the sigmoidal tone 
curve 1 70 asymptotes to 0 and 1 00 at negative and positive infinity. Therefore, once the tone cun/e 1 70 is generated, 
the curve 170 must be scaled to th6[0, 100] range after to make sure 0 maps to 0 and 100 maps to 100 as follows: 

M^= ^SVff - 

[0032] The sigmoidal shaped tone curve 170 of the type shown in FIGURE 9 that is produced by the tone curve 
generator 160 has several properties that satisfy the constraints and perceptual preferences discussed above. The 
sigmoidal shaped tone curve 170 is mostly linear at the mid range, with saturation at the dark and bright end if the 
slope Is high. A mostly linear segment at the mid range means that for most Images, the bulk of the pixels will be 
mapped more or less linearly, thus presen/ing the shape of the original histogram. In addition, the parameters of the 
sigmoidal tone curve 170 can be set to achieve the desired standard deviation level. Furthermore, the sigmoidal 
function is smooth, with rather consistent curvature across different parameter settings (accelerating before the p point, 
decelerating after). Therefore, the sigmoidal function is unlikely to generate tone curves 1 70 that have wild bends and 
50 peaks. 

[0033] Once the tone curve 170 is generated (step 240), the received linear image data 130 is mapped onto this 
tone curve 170 by the tone mapping operator 40 to produce output image data 45 corresponding to the levels of the 
output device 50 (step 250). It should be noted that the output device 50 can be a display, printer or other image 
processing device, it should also be noted that additional processing may be performed on the output image data 45 
^ prior to being sent to the output device 50 (step 260). 

[0034] With reference now to FIGURES 3A-3C of the drawings and the corresponding flow charts of FIGURES 4A- 
4Cj alternative embodiments for implementing the image statistics processor of FIGURE 1 are illustrated. In a first 
embodiment, as shown in FIGURE 3A and described in FIGURE 4A, the linear image data 130 (step 400) is passed 



5 



EP 1 302 898 A2 



to a histogram generator 122, which produces a histogram 145 (e.g.. a table) describing the distribution of the digital 
values of the data 130 (step 405). For example, In a color digital Image, there are a certain number of red, blue and 
green values (e.g., one of each for each pixel in the image). In order to produce a histogram of green values, for 
example, the histogram generator 122 determines a bin size (e.g.. range of values, such as 0-10, 11-20, 21-30, etc.) 
and counts the number of occurrences of green pixel values in each bin, To create a graphical representation of the 
histogram, the data In each bin Is plotted as a function of the bin centers (e.g., 5, 15, 25, etc.). 
[0035] Once the histogram is created, the histogram bin center values 135 are transmitted to an L * bin center 
converter 1 24 to transform the linear bin center values 1 35 into L * bin center values 1 35a (step 410) The conversion 
is performed using equation (2) above, with y being the bin center values, and being the largest possible digital 
value for the pixel values, The L * bin center values 136a and the histogram 145, which Includes the count of the 
number of values in each bin are transmitted to L* statistics calculation logic 128, which calculates the Image statistics 
149 (step 420) that are subsequently used in the calculation of the slope parameter cc and the shift parameter p for the 
sigmoidal tone curve. 

[0036] The various embodiments for calculating the slope parameter ct and the shift parameter p are discussed below 
with reference to FIGURES 5 and 6. Briefly, the Z_* image statistics 149 used for calculation of the slope parameter a 
and the shift parameter p include the current L * standard deviation and either the current mean L * value or the L * 
percentile values of the histogram (Mh and (1 00-n)th percentiles). The Z_* percentile values of the histogram represent 
the end points of the histogram. Therefore, n Is chosen to be relatively small (preferably less than 5). For example, If 
5% of all L * pixel values have values less than 18, then the 5^^ percentile is 18, thus the low end point would be 1 8 . 
The high end point would then Just be the L * value corresponding to the 95*"^ percentile (e.g., 95% of ail L * pixel 
values have values lower than the high end point), 

[0037] The current L * standard deviation o^^j^f^fOeLn be calculated from the histogram as follows: 



where BfS are histogram bin centers In L * values, /?/ Is the pixel count In the Ith bin of the histogram, n^, Is the number 
of histogram bins, n is the total number of pixels In the Image, and I Is the mean pixel L « value, which can also be 
calculated from the image histogram as follows: 



[0038] In one embodiment, c^current 's calculated from the histogram on the green color plane only. For most sensors, 
the green channel Is the color plane that approximately captures the luminance Information. However, It should be 
noted that some combination of the RGB pixel values can be used instead to calculate the image L* standard deviation. 
The resulting tone curve should not vary significantly as long as the choice Is a reasonable representation of the 
luminance. 

[0039] Alternatively, as shown in FIGURE 3B and described in FIGURE 4B, if the histogram is not available, the 
Image statistics 1 49 can be calculated from the Image pixel values themselves. Therefore, when linear Image data 1 30 
arrives at the image statistics processor 1 20 (step 430), the linear image data 1 30 is converted into L * image values 
1 30a by an L * image data converter 1 25 (step 435), and the L * image values 1 30a are passed onto the L * statistics 
calculation logic 128, where the image statistic 149 (e.g., current L * standard deviation and either the current mean 
L * value orthe L ^ percentiles) are calculated (step 440). For example, in this embodiment, the current mean L * value 
can be calculated by taking the average of all of the L image values. 

[0040] As a further alternative, as shown in FIGURE 3C and described in FIGURE 4C, to reduce computational 
complexity, since the histogram bins do not change from Image to Image (only the courts In each bin change), the 
histogram bin center values 1 35 can be converted into L * bin center values 1 35a prior to receiving the linear image 
data 130. Therefore, in this embodiment, once the histogram bin center values 135 are determined (step 450), the bin 
center values 135 are converted to L * bin center values 135a by the L * bin center converter 124 (step 455) and these 
L * bin center values 1 35a are stored in a memory 129 accessible to or within the image statistics processor 1 20 (the 
latter being illustrated] (step 460). As a further alternative, theconversion to L>f bin center values 135a can be performed 
externally, and the L * bin center values 135a can be pre-loaded into the memory 129. 




(Equation 4) 




(Equation 5) 
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[0041 ] Thereafter, when the linear Image data 1 30 arrives at the Image statistics processor 1 20 (step 465), the linear 
Image data 1 30 is sent to the histogram generator 1 22 (described above), which counts the number of pixel values In 
each bin (step 470). Once counted, the table 145, which Includes the count of the number of pixel values in each bin 
is sent to the L * statistics caicuiation logic 12B, which retrieves the L *bin center values 135a from the memory 129 

5 (step 480) to calculate the image statistics 149 (step 485). With this embodiment, the image statistics 149 {e.g., L * 
standard deviation and mean L* value) can be calculated with only simple addition, subtraction and multiplication. 
[0042] With reference now to FIGURES 5, 6A and 6B, there is illustrated the steps for calculating the slope parameter 
a and the shift parameter p using the image statistics (e.g., L * standard deviation and either mean L * value or L 
^percentiles) calculated by the image statistics processor. In calculating the shift parameter and the slope parameter, 

^0 as discussed above, various perceptual preference guidelines can be used to produce an output image that is pleasing 
to the viewer. The perceptual preference guidelines can be chosen by the user or can be selected by the manufacturer 
of the specific digital image system in which the tone mapping algorithm is implemented. As an example, the following 
discussion is based on three preferred perceptual preferences, including the presen/ation of general histogram shape 
at the middle range, a desired L * standard deviation level of around 20 and a desired mean L * level of around 50. It 

15 should be understood that other perceptual preferences can be used instead of these preferred preferences. 

[0043] As illustrated In FIGURE 5, to determine the slope parameter oc, the current L * standard deviation calculated 
by the image statistics processor is transmitted to the tone curve parameter caicuiation logic (shown in FIGURE 1) 
(step 600). The desired L * standard deviation is pre-selected by the user or the manufacturer and stored within a 
memory (not shown) accessible to or within the tone curve parameter calculation logic (step 51 0). As discussed above, 

20 in preferred embodiments, the desired L * standard deviation is chosen to be about 20. However, other values can be 
chosen, depending upon the requirements of the user or the manufacturer Values greater than 20 produce a more 
contrasty image, while values less than 20 produce a less contrasty image, in addition to the desired L standard 
deviation, a default slope value Is also pre-selected (step 520). 

[0044] With a linear tone curve, the standard deviation can be changed from one value (p) to another value [q) by 
25 using a tone curve slope of q/p. The sigmoidal tone curve is nonlinear, but the mid-section of the sigmoidal curve is 
mostly linear Since for most images, a large proportion of image pixels have values close to the mid-range (around 
the point x = L for one of the embodiments), the slope of the sigmoidal curve can be changed at this point to try to 
achieve the desired L * standard deviation. 

[0045] At the point X = 1 0Op, the slope of t(x} Is 1 when a = 4. Therefore, the default slope a,^^^Jfi can be set to 4. 
50 However, the actual tone curve used (as discussed above) Is t^{x), which has a slope of 1 when a is less than 4. In 

practice, any value between 3 and 4 can be used as the default slope ^^efaufb depending on how contrasty one would 
iii<e the image to look. Once the default slope has been selected, the slope parameter cc can be calculated based on 
the current L standard deviation, the desired £.* standard deviation and the default slope (step 530). To get the slope 
parameter cc for an Image with an L * standard deviation of ctcunsnt' default a value is scaled according to the current 
35 standard deviation and the desired L * standard deviation o^.-,„w as follows: 



a = a^efauft X (Equation 6) 

^ current 

40 

[0046] As mentioned above, the shift parameter p shifts the histogram to affect the general brightness of the output 
image. There are different ways to determine the shift parameter, depending on what perceptual factors are more 
important in the particular application, and what computational complexity Is allowed. FIGURES 6A and 6B Illustrate 
two such alternative embodiments for calculating the shift parameter p. 

[0047] In a first embodiment, as shown In FIGURE 6A, the shift parameter |3 Is calculated based on the mean L >« 
value. In this embodiment, the current mean L * value calculated by the Image statistics processor Is sent to the tone 
curve parameter calculation logic (shown in FIGURE 1) (step 600). The desired mean value Is pre-selected by the 
user or the manufacturer and stored within a memory (not shown) accessible to or within the tone cun/e parameter 
caicuiation logic (step 610). Under some circumstances, however, this pre-selected desired mean L * value may be 
50 adjusted based on additional Information about the image being processed. For example, If the content of the image 
is known to the image processor to be predominantly black objects, then the desired L * value can be adjusted to be 
lower than 50, to avoid over-brighten In g the Image. Or, If the Image data Is known to have a high level of noise, the 
desired L ^ value can be biased to be closer to the original Image mean L ^ value, to avoid amplification of the noise 
through tone mapping. 

55 [0048] For example, the desired mean L * value of the whole image can be pre-selected to be around 50 (medium 
brightness). Alternatively, the desired mean L * value can be selected to be 50 in a certain area of the Image. As an 
example, for images containing a human face, the L ^ mean value at the face area can be set to 50. Using such criterion 
to determine the shift parameter of a sigmoidal tone curve generally results In a good quality image. However, when 
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face recognition Is not feasible, the nnean L * value of the center part of an Image, or some other area of Interest In the 
Image, can be brought to a level of around 50. This latter method of determining the shift parameter Is more practical 
for a low cost Imaging system. It should be understood that any other desired mean L * value can be used instead of 50 . 

[0049] The shift parameter can be calculated based on the current mean L * value and the desired mean L * value 
5 (step 620). For example, if L is the mean L ^ value of the area of interest in an image (e.g., a face area), and L^^simd 
Is the desired mean L * value, the shift parameter Is equal to: 

P = (L - L^i^)nQQ + 0.5. (Equation 7) 

10 

When using an L^^g^,-^^^ value of 50, p is simply L/100. 

[0050] In a second embodiment, shown in FIGURE 6B, the shift parameter p is calculated based on the L ^ percentile 
values. Various studies ("se^, e g., Tanaka et al., 1997) have suggested that obsen/ers preferred renderings that have 
as large an output gamut as possible. Therefore, to accommodate these preferences, the shift parameter can be 

t5 determined by shifting the histogram to mal<e it occupy as much of the output device dynamic range as possible, while 
bringing as little disturbance to the shape of the histogram as possible to avoid saturation of a lot of points at either 
end of the distribution. This is done by using a centering function in which the tone curve is centered at midway between 
the end points of the histogram. This Is a conservative choice, typically resulting In only small changes to the mean 
luminance of the image if the original RGB histogram is somewhat "normal" (values are not concentrated on a small 

20 section of thefuil range). 

[0051] Therefore, when the tone curve parameter calculation logic receives the percentiles (e.g., nth and (100 - 
n)Vr\ percentiles, with n being a convenient small value) (step 650), the shift parameter |3 can be calculated as the mid- 
point between these two end points (step 660). For example, If the nth percentile of pixel values In the histogram Is x.|, 
and the (100-n}th percentile of pixel values Is Xg, then the shift parameter is determined by the following equation: 

25 

^ ^ + Xg (Equation 8) 

SO The combined histogram of RGB values Is preferably used to determine the percentiles, so that no color plane Is 

ignored when calculating the shift of the tone curve, This is especially Important for tone mapping performed on Images 
with a strong color cast, to avoid color shifts after tone mapping. 

[0052] Both methods of determining p discussed above have been implemented and tested on Images. The first 
method (P calculated from mean L * value) is simpler to implement, and works well on most Indoor and outdoor Images. 
35 When the mean L * value of the whole image is used as p. It tends to over-brighten night-time or dark images, i.e., 
images that are meant to have a low mean value, such as candlelight scenes. The second method (using the average 
of the top and bottom percentile of the histogram as p) works better for dark images. It tends to change the image 
brightness less than the first method^ so it is a more conservative method. 

[0053] For some applications, it may be desirable to iimitthe extent that the sigmoidai curve is shifted, i.e., the extent 
40 the shift parameter p can deviate from 0.5. For example, for severely under-exposed Images, excessive shifting up of 
the tone curve (P value close to 0) can result in a low-quality image due to amplification of sensor noise. Likewise, 
excessive shifting down of the tone cun/e (P value close to 1 ) for highly over-exposed Images can produce color shifts 
in nearly saturated areas. To minimize these problems, the value of the shift parameter can be limited depending on 
the mean pixel values of the original image, if the original image mean is too high or too low, depending on the specific 
application of the tone mapping algorithm, then the shift parameter p Is adjusted to deviate less from the midpoint 0.5 
In order to reduce noise amplification. In one embodiment, a series of p value upper limits 7} are pre-selected, each 
corresponding to one mean level Mj of the original image pixel values. When the original image mean pixel value is 
less than Mp then the p value is limited to be at least Tj, i.e., if the value calculated from Equation 7 or Equation 8 is 
less than 7], then p is set to be equal to Ty. Similarly, a series of lower limits 7"^ can also be pre-selected^ each corre- 
ct? spending to one mean pixel level M^. When the original image mean pixel value is larger than M^^ then p Is limited to 
be no greater than By limiting the shift parameter in this way, artifacts resulting from tone over-compensation can 
be minimized. 

[0054] With reference now to FIGURES 7A and 7B, there are illustrated alternative embodiments for implementing 
the tone curve generator 160 shown in FIGURE 1 of the drawings. In a first embodiment, as shown In FIGURE 7A and 
55 described in the steps of FIGURE 8A, when the tone curve parameters 150 (e.g., slope parameter a and the shift 
parameter p) have been calculated by the tone curve parameter calculation logic 140 (shown in FIGURE 1), these 
parameters are sent to the tone curve generator 160 (step 800). Thereafter, the sigmoidai tone curve 1 70 is generated 
by sigmoidai function logic 1 62 within the tone curve generator 1 60 using the received tone curve parameters 1 50 and 
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equations (1) and (3) listed above (steps 805 and 810), Once generated, the resulting tone curve 170 resembles the 
tone curve shown In FIGURE 9. 

[0055] The tone curve 170 shown in FIGURE 9 is capable of mapping L * image input values onto L* image output 
values forthe appropriate output device. However, since L * is not acommonly used scale n output devices, the output 
5 digital image data is converted to linear form before conversion to the scale required by specific output devices. In 
preferred embodiments, the Y-axis (L * Image output values) Is converted by a Y-axIs converter 1 64 from an L * scale 
to a linear scale (step 815), which is approximately a power of three transformation, as discussed above, to produce 
tone curve 1 70a. 

[0056] In addition, since most display devices, such as a CRT display, have an Intensity to voltage response curve 

^0 that is roughly a power function (e.g., If the pixel value is set at x, the displayed luminance Is proportional to with y 
being the display's gamma value), the linear output values need to be gamma corrected in orderto be displayed properly 
on such output devices. 1 n most tone mapping methods; gamma correction is performed after mapping. In the sigmoidal 
tone mapping method described herein, gamma correction can either be performed separately after tone mapping, or 
can be combined with the tone mapping step. In one embodiment, the gamma correction is performed together with 

IS tone mapping by applying gamma correction logic 1 66 to the Y-axis to produce tone curve 1 70b (step 820). The com- 
bined effect of linear transformation and gamma correction of tine Y-axis (output) is a power of 3/7, where y is the 
intended display gamma value, such as 2.4, if the output device Is an sRGB (standard color space) display. 
[0057] Typically, the tone mapping Is applied to linear RGB Input pixel values (e.g., linear image data 130 Is supplied 
by the linear transformation device 30 to the tone mapping operator 40 in FIGURE 1), The sigmoidal tone curve is 

20 generated assuming an Z_ * scale on the X-axis, I.e., X-axis values uniformly sampled on the L scale. Therefore, In 
order to properly map linear input values onto the tone cun/e, the tone curve's X-axIs sampling has to be converted to 
be uniform on a linear scale for easy table iool<up operation when applying the tone cun/e to images. This is done by 
an X-axIs converter 168 to produce tone curve 170c (step 825), using a simple linear Interpolation process, which 
requires only arithmetic operations, as described above In connection with the Y-axis converter. This final tone curve 

25 170c is used by the tone mapping operator (shown in FIGURE 1) to map the linear image data onto the luminance of 
the output device. 

[0058] In a second embodiment, as shown in FIGURE 7B and described in the steps of FIGURE SB, since one of 
the most computationally costly parts of the tone mapping algorithm Involves the exponential function calculation in 
the sigmoidal function, to reduce computational complexity, an approximation of the sigmoidal tone curve calculation 

50 can be Implemented with only arithmetic operations. The exponential calculation In the sigmoidal function can be avoid- 
ed by pre-calculating a pair of fixed tone cun/es with two extreme slopes, and then interpolating to obtain tone curves 
with Intermediate slopes. As discussed above, the shift parameter translates the tone curve on the horizontal (X) axis. 
Therefore, by pre-calculating the tone cun/es for a wider range on the X (input pixel value) axis, appropriate tone curves 
can be selected by "cutting out" different sections of the pre-calculated tone curves for different shift parameters. 

35 [0059] To pre-calculate the two tone curves, the user or manufacturer must pre-select two slope parameters, a^^x 
and ttf^jn (step 850). The slope parameter a for most images falls in a range from 2.5 to 6.5. Therefore, in preferred 
embodiments, ctmin set to be 2.5 and a^^x Once the maximum and minimum slope parameters 1 55 

are set, sigmoidal function logic 162 can pre-calculate the two tone cun/es (step 855), one with a slope of 2.5 (curve 
P^), which gives an almost linear tone curve (on the L * scale), and one with a steep slope of 6.5 (cun/e P2), which 

40 gives a lot of contrast enhancement to the Image. The two pre-generated tone curves are shown In FIGURE 10. It 
should be noted that the two curves are generated for a shift value p of 0, (The X-axis range for different p values are 
later determined by translating the X origin on the X-axis, as discussed hereinbelow.) 

[0060] For slope values between 2.5 and 6.5, the tone curve is generated by linear interpolation between and Pg. 
Let Pa(x) be the tone curve with an a value of a, with a c [2.5, 6.5], then : 

P^(x) = P,{x) + (P2(x) - P,(x)) *g^^. (Equation 9) 

50 [0061 ] This of course is an approximation, but for the limited a value range between 2.5 and 6.5, the approximation 
is reasonably accurate, in FIGURE 11 A, the tone curve output L * values are plotted as a function of a for 11 different 
input pixel levels (represented by the 11 lines). As can be seen^ the output L * values are close to linearly related to a. 
[0062] As discussed above, typically, the tone mapping output is gamma-corrected linear RGB values ready for 
display on a screen. Therefore, as described above in connection with FIGURES 7A and 8A, the Y-axis of P^ and P2 

55 158 may be converted to gamma-corrected linear RGB values (steps 860 and 865) by the Y-axis converter 1 64 and 
gamma correction logic 166 of the tone curve generator 160. When the tone curve output Pg[x) is gamma corrected 
linear RGB values, the relationship with oc is still quite linear, as shown In FIGURE 11 B. Therefore, by pre-calculating 
P^ and P2 158 as gamma corrected pixel values, the Interpolated tone values will be gamma corrected as well. This 
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means the tone curve calculation only requires the arithmetic operations used In Interpolation (equation 9 above). Once 
the gamma corrected tone curves Pi and P2 158 have been pre-generated, they are stored In a memory 169 (step 
870) accessible to or within the tone curve generator for later use in calculating the tone curve for a particular set of 

input values. 

5 [0063] When the tone curve generator 1 60 receives the current slope parameter CLj^urrent current shift 

parameter Pcurrent "^^^^ (st^P 875) from the tone curve parameter calculation logic (shown In FIGURE 1), calculation 
logic 163 within the tone curve generator 160 determines Pa(^ "I^S using the gamma corrected tone curves P-, and 
P2 158 and the current slope parameter 150a and equation 9 above (step 880). Thereafter, to select a section of 
(X), the current shift parameter 150b and Pa(x) 159 are transmitted to shifting logic 165 within the tone curve generator 

^0 1 so (step 885). Since the shift parameter 150b simply translates the tone curve on theX-axIs, tone curves with different 
shift parameters can be obtained by translation on the X-axis of Pq(x). 

[0064] For example, assuming PJ^x) is generated with a shift parameter of 0 and in the range [-1 00, 100], to get a 
tone curve with a shift parameter of |3 at the X-axis range of [0, 100], the tone curve is translated on the X-axis so that 
an X value of - p 1 00 now becomes 0, and an X value of (1 - P) * 1 00 becomes 1 00. The tone curve with shift parameter 
IS p is obtained by "cutting out" the pre-generated tone curve Pg(x)^rom the point X= - p * 100 to the point X= (1 - P) * 
1 00, In Impiementation, this operation involves simply moving the file pointer to different starting points in a pre-gen- 
erated table, thus requiring almost no computation. 

[0065] After the tone curve Is "cut out" from P^ (x), the two ends of the tone curve need to be scaled to the full Y 

range. This again only requires arithmetic operations on each entry of the tone mapping table that was previousiy cut 
20 out from the pre-generated table. Thereafter, as discussed above In connection with FIGURE 7A and 8A, the X-axis 
converter 1 68 converts the X-axIs from an L --^ scale to a linear scale (step 890) to produce the final tone curve 170. it 
should be noted that the translation of the tone curve can happen before or after the interpolation step required to get 
the correct slope by translating tone curves P^ and Pg 158. For the least amount of computation, the translation Is 
performed first, so that the Interpolation only needs to be done on the "cut out" section of the pre-generated tone curves 
25 and P2 158. 

[0066] As wili be recognized by those skilled In the art, the Innovative concepts described In the present application 
can be modified and varied over a wide range of applications. Accordingly, the scope of patented subject matter should 
not be limited to any of the specific exemplary teachings discussed, but is instead defined by the following claims. 

so 

Claims 

1. A digital Image system (100) capable of receiving digital Image data (25) representing at least part of a digital 
Image and mapping said digital Image data (25) onto the range of an output device (50), said digital Image system 

35 (100) comprising: 

an image statistics processor (120) capable of calculating at least one current image statistic (149) based on 
said digital image data (25); 

calculation logic (140) connected to receive said at least one current image statistic (149) from said image 
40 statistics processor (120) and further being capable of calculating at least one tone curve parameter (150) 

based on said at least one current image, statistic (149) and at least one perceptual preference associated 

with said digital Image system (100); and 

a tone curve generator (160) connected to receive said at least one tone curve parameter (150) from said 
calculation logic (1 40) and further being capable of generating a tone curve (1 70) using said at least one tone 
curve parameter (150) and a sigmoldai function, said tone curve (1 70) being used to map said digital Image 
data (25) onto the range of said output device (50). 

2. The system (1 00) of Claim 1 , further comprising: 

50 a linear transformation device (30) connected to receive said digital image data (25), convert said digital image 

data (25) into linear image data (130) and transmit said linear image data (130) to said image statistics proc- 
essor (120), said Image statistics processor (120) using said linear Image data (130) In calculating said at 
least one current Image statistic (1 49). 

55 3. The system (100) of Claim 1 or 2, wherein said image statistics processor comprises: 

a histogram generator (122) connected to receive said linear image data (130) and generate a histogram (145) 
of the pixel values represented by said linear image data (135), said histogram (145) containing a plurality of 
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bins representing respective ranges of pixel vaiues and a count of the nunnber of pixei vaiues represented by 
said ilnear image data (130} witliin eacii of said piuraiity of bins; 

a converter (124) capabie of converting a bin center vaiue (135) for eacli of said bins associated witli said 
histogram (1 45) from a linear scale to an L * scale to produce a plurality of L * bin center values (1 35a); and 
5 statistics calculation logic (128) connected to receive said count from said tiistogram generator (122) and said 

piuraiity of L * bin center vaiues (1 35a) from said converter (1 24) and caicuiate said at ieast one current image 
statistic (149) using said count and said piuraiity of L bin center vaiues (135a). 



4. Tlie system (1 00) of Ciaim 3, furtiier comprising: 

10 

a memory (129) for storing said piuraiity of L * bin center vaiues (135a) prior to said iiistogram (145) being 
generated, said statistics calculation iogic (128) retrieving said piuraiity of L » bin center vaiues (135a) in 
response to receiving said count. 

t5 5. The system (100) of one of the preceding Claims, wherein said at least one tone curve parameter (1 50) Includes 
a slope parameter and a shift parameter and wherein said at least one current Image statistic (149) includes a 
current L* standard deviation and one of a current mean L * vaiue or iow and high L * percentile vaiues. 

6. The system (100) of Ciaim 5, wherein said at least one perceptual preference includes a desired L standard 
20 deviation and one of a desired mean L * value or a centering function, said slope parameter being calculated using 

said current L * standard deviation and said desired L standard deviation, said shift parameter being calculated 
using either said current mean L ^ vaiue and said desired mean L * value or said iow and high L percentile values 
and said centering function. 



25 7, The system (100) of Claim 6, further comprising: 



an upper pre-selected mean L * value and an associated upper pre-seiected shift vaiue, said shift parameter 
being set to said upper pre-seiected shift vaiue when said current mean L * vaiue is less than said upper pre- 
seiected mean vaiue and said caicuiated shift parameter is less than said upper pre-seiected shift value; and 
so a iower pre-seiected mean L * vaiue and an associated iower pre-selected shift vaiue^ said shift parameter 

being set to said lower pre-seiected shift vaiue when said current mean L ^ value is greater than said lower 
pre-seiected mean L * vaiue and said caicuiated shift parameter is greater than said iower pre-selected shift 
vaiue. 



55 8. The system (1 00) of one of Claims 5 - 7, further comprising a memory (169) for storing first and second pre- 

calcuiated tone curves (1 58) generated by said tone curve generator (1 60) prior to said digital image system (1 00) 
receiving said digital image data (25), said first p re-calculated tone curve having a minimum slope and said second 
pre-caiculated tone curve having a maximum siope. 



40 9. The system of Ciaim 8, wherein said tone curve generator comprises: 



calculation logic (1 63) connected to receive said siope parameter (1 50a) and said first and second pre-calcu- 
lated tone curves (158), said calculation logic (163) being further capable of interpolating between said first 
and second pre-caiculated tone curves (158) using said siope parameter (150a) to obtain an initial tone curve 
^ (159); and 

shifting iogic (165) connected to receive said initiai tone curve (159) and said shift parameter (150b), said 
shifting iogic (1 65) being further capable of shifting said initial tone curve (1 59) on the X-axis using said shift 
parameter (150b) to produce said tone curve (170) used in mapping said digitai image data (25) onto the range 

of said output device (50). 

50 

10. The system of one of the preceding Ciaims, wherein said tone cun/e generator further comprises: 



gamma correction logic (1 66) for applying gamma correction to the Y-axis of said tone curve (170); and 

a converter (1 68) for converting theX-axIs of said tone curve (1 70) to the scale of said digital image data (25) . 

55 

11. A method for mapping digital image data (25) representing at least part of a digitai image onto the range of an 
output device (50), said method comprising: 
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receiving said digitai image data (25} at a digltai image system (100); 
caiculating at ieast one current image statistic (149) based on said digltai image data (25); 
calculating at least one tone curve parameter (150) based on said at ieast one current image statistic (149) 
and at least one perceptual preference; and 
5 generating a tone curve (170) using said at least one tone curve parameter (150) and a sigmoldal function, 

and using said tone curve (1 70) to map said digitai image data (25) onto tlie range of said output device (50). 

12. Tine method of Claim 11 , wherein said step of calculating said at least one current image statistic (149) further 
comprises: 

10 

converting said digital Image data (26) Into linear image data (130); and 

caiculating said at ieast one current image statistic (149) using said linear image data (130), 

13. The method of Claim 12, wherein calculating said at least one current image statistic (149) further comprises: 

15 

generating a histogram (145) of the pixel values represented by said linear Image data (130), said histogram 
(145) containing a plurality of bins representing respective ranges of pixel values and a count of the number 
of pixel values represented by said linear image data (130} within each of said plurality of bins; 
converting a bin center value (135) for each of said bins associated with said histogram (145) from a linear 
20 scale to an L * scale to produce a plurality of L bin center values (135a); and 

calculating said at ieast one current image statistic (149) using said count and said plurality o\ L* bin center 
values (135a). 

14. The method of Claim 13, wherein converting said bin centers further comprises: 

25 

converting said bin center values (135) to produce a plurality of L « bin center values (1 35a) prior to said step 
of generating; and 

storing said plurality of /. * bin center values (135a) within a memory (129). 

so 15. The method of one of Claims 11 - 14, wherein said at least one tone curve parameter (150) includes a slope 

parameter and a shift parameter and wherein said at ieast one current image statistic (149) includes a current L 
* standard deviation and one of a current mean L * value or low and high L * percentile values. 

16. The method of Claim 15, wherein said at least one perceptual preference includes a desired L * standard deviation 
35 and one of a desired mean L value or a centering function, calculating said at least one tone curve parameter 

(150) further comprising: 

caiculating said slope parameter using said current L « standard deviation and said desired L ^ standard 
deviation; and 

calculating said shift parameter using either said current mean L value and said desired mean L * value or 
said low and high L f percentile values and said centering function. 

17. The method of Claim 16, wherein calculating said at least one tone curve parameter (150) further comprises: 

setting said shift parameter to be equal to an upper p re-selected shift value when said current mean L « value 
is less than an upper pre-selected mean L ^ value associated with said upper pre-selected shift value and 
said calculated shift parameter is less than said upper pre-selected shift value; and 

setting said siilft parameter to be equal to a lower pre-selected shift value when said current mean L * value 
is greater than a lower pre-selected mean L ^ value associated with said lower pre-selected shift value and 
50 said calculated shift parameter is greater than said lower pre-selected shift value. 

18. The method of one of Claims 15-17, wherein generating further comprises: 

generating first and second pre-caiculated tone curves (158) prior to said step of receiving said digltai image 
55 data (25), said first pre-caiculated tone curve having a minimum slope and said second pre-caiculated tone 

curve having a maximum slope. 

19. The method of Claim 18, wherein generating said tone curve (170) used in mapping said digital image data (25) 
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onto the range of said output device (50) comprises: 

interpolating between said first and second p re-ca leu late d tone curves (158) using said slope parameter (1 50a) 
to obtain an initial tone curve (159); and 
5 shilting said initial tone curve (1 59) on the X-axis using said shift parameter (1 50b) to produce said tone curve 

(170) used in mapping said digital Image data (25) onto the range of said output device (50). 

20. The method of one of Claims 11 - 19, wherein generating said tone curve (170) comprises: 

^0 applying gamma correction to the Y-axis of said tone curve (170); and 

converting the X-axIs of said tone curve (1 70) to the scale of said digital Image data (25). 

15 
20 
25 

so 

35 
40 
45 
50 
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